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Abstract 


The generalized mathematical formulation of sloshing dynamics for partially 
filled liquid of cryogenic superfluid helium II in devar containers driven by 
both the gravity gradient and jitter accelerations applicable to scientific 
spacecraft which is eligible to carry out spinning motion and/or slew motion for 
the purpose to perform scientific observation during the normal spacecraft 
operation are investigated. An example is given with Gravity Probe-B (GP-B) 
spacecraft which is responsible for the sloshing dynamics. The jitter 
accelerations include slew motion, spinning motion, atmospheric drag on the 
spacecraft, spacecraft attitude motions arising from machinery vibrations, 
thruster firing, pointing control of spacecraft, crew motion, etc. Explicit 
mathematical expressions to cover these forces acting on the spacecraft fluid 
systems are derived. The numerical computation of sloshing dynamics have been 
based on the non- inertia frame spacecraft bound coordinate, and solve time- 
dependent, three-dimensional fprmulations of partial differential equations 
subject to initial and boundary conditions. The explicit mathematical 
expressions of boundary conditions to cover capillary force effect on the liquid- 
vapor interface in microgravity environments are also derived. The formulations 
°f fluid moment and angular moment fluctuations in fluid profiles induced by the 
sloshing dynamics, together with fluid stress and moment fluctuations exerted on 
the spacecraft dewar containers have been derived. Results of this study have 
been widely published in the open journals. 
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I. Introduction 


For the purpose to carry out scientific experiments, some experimental 
spacecraft use cryogenic cooling for observation instrumentation and telescope, 
superconducting sensors for gyro read-out and maintain very low temperature near 
absolute zero for mechanical stability. The approaches to both cooling and 
control involve the use of superfluid liquid helium II. In this study, sloshing 
dynamics associated with spinning motions are investigated. To cover the 
spacecraft spinning motions, the Gravity Probe-B (GP-B) spacecraft has been 
chosen as the example in this study. The GP-B spacecraft adopts the cooling and 
boil -off from the cryogenic liquid helium dewar as a cryogen and propellant to 
maintain the cooling of instrumentations, attitude control and drag-free 
operation of the spacecraft. The potential problems for cryogenic liquid in 
dewar container could be due to asymmetry in the static liquid helium 
distribution and to perturbations in the liquid-vapor interface caused by slosh 
wave excitation driven by pointing control, machinery vibration, etc. 

For the case of the GP-B spacecraft, cryogenic liquid helium II, at a 
temperature of 1.8 K, is used as the propellant. With its superconducting 
behavior, there is no temperature gradients in the liquid helium. In the absence 
of temperature gradient along the surface which drive Marangoni convection, the 
equilibrium shape of the free surface is governed by a balance of capillary, 
centrifugal and gravitational forces. Determination of liquid-vapor interface 
profiles based on computational experiments can uncover details of the flow which 
can not be easily visualized or measured experimentally in a microgravity 
environment . 

The instability of the liquid-vapor interface surface can be induced by the 
presence of longitudinal and lateral accelerations. Slosh waves are, thus, 
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excited which produces high and low frequency oscillations in the liquid 
propellant. The sources of the residual accelerations range from the effects of 
the Earth's gravity gradient and jitter accelerations which include, atmospheric 
drag on the spacecraft, vibration of compressor, spacecraft attitude motions 
arising from machinery vibrations, thruster firings, spacecraft slew motion, 
pointing control of spacecraft, crew motion, etc. Recent study (Kamotani et 
al. , 1981) suggest that the high frequency accelerations may be unimportant in 
comparison to the residual motions caused by low frequency accelerations. 

Time - dependent dynamical behavior of partially- filled rotating fluids in 
reduced gravity environments was simulated by numerically solving the Navier 
Stokes equations subject to the initial and the boundary conditions (Hung and 
Shyu, 1991 a,b,c; 1992 a,b,c; Hung et al., 1991 a,b,c; 1992 a,b,c). At the 
interface between the liquid and the gaseous fluids, both the kinematic surface 
boundary condition, and the interface stress conditions for components tangential 
and normal to the interface, were applied (Hung and Shyu, 1991 a,b,c; 1992 a,b,c; 
Hung et al., 1991 a,b,c; 1992 a,b,c). The initial conditions were adopted from 
the steady-state formulations developed by Hung et al (1989 a,b,c). Some of the 
steady- state formulations of interface shapes were compared with the available 
experiments carried out by Leslie (1985) in a free-falling aircraft (KC-135). 
The experiments carried out by Mason et al (1978) showed that the classical fluid 
mechanics theory is applicable for cryogenic liquid helium in large containers. 

In the spacecraft orbit around the Earth, the direction of azimuth angle 
of Earth toward the location of the spacecraft geometric center varies from 0° 
along the rolling axis of spacecraft to various directions in which three 
dimensional calculation shall be assumed. 

As the spacecraft moves along the orbit, any fluid capable of motion 
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relative to the spacecraft is subject to the acceleration that arises from the 
gravity gradients of the Earth (Avduyevsky, 1984; Forward, 1982; Misner et al. , 
1973) . The interaction between the particle mass of fluid and the spacecraft 
mass due to gravity gradient accelerations (Forward, 1982) are capable for the 
excitation of slosh waves and disturb the fluid system which induces the 
fluctuations of viscous stress and its moment exerted on the containers of the 
spacecraft fluid systems. In the meanwhile, the sources of residual acceleration 
of gravity jitter range from atmospheric drag on the spacecraft, background 
gravity, spacecraft attitude motions arising from machinery vibrations, 
spacecraft slew motion, thruster firings, crew motion, etc. , are also capable for 
the excitation of slosh waves on the fluid containers. 

It is critically important to understand the physical and dynamical 
behavior of cryogenic helium in a rotating cylinder to effectively promote space - 
oriented missions. 

At temperatures close to absolute zero, quantum effects begin to be of 
importance in the properties of fluids. At a temperature of 2.17°K, liquid 
helium has a A-point (a second-order phase transition); at temperatures below 
this point liquid helium (helium II) has a number of remarkable properties, the 
most important of which is superfluidity. This is the property of being able to 
flow without viscosity in narrow capillaries or gaps. 

The basis of the dynamics of helium II is the following fundamental result 
of microscopic theory. At temperatures other than zero, helium II behaves as if 
it were a mixture of two different liquids . One of these is a superfluid and 
moves with zero viscosity along a solid surface. The other is a normal viscous 
fluid. It is of great importance that no friction occurs between these two parts 
of the liquid in their relative motion, i.e. , no momentum is transferred from one 
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to the other. 


It should be emphasized that regarding the liquid as a mixture of normal 
and superfluid parts is no more than a convenient description of the phenomena 
which occur in a fluid where quantum effects are important. One of these motions 
is normal and has the same properties as the motion of an ordinary viscous fluid, 
but the other is the motion of a superfluid. The two motions occur without any 
transfer of momentum from one to another. We can, in a certain sense, speak of 
the superfluid and normal parts of the fluid, but this does not mean that the 
fluid can actually be separated into two such parts. 

If fluid flow can separate helium II into the regions of the superfluid and 
the normal fluid, two temperatures zones are immediately created. A very low 
temperature zone is located at the zone of very high density concentration of the 
superfluid, while a high temperature (below 2.17° K) zone is located at the zone 
of very high density concentration of the normal fluid at the ocher end. The 
existence of a sharp temperature gradient at the interface between the superfluid 
and the normal fluid results in the creation of a great difference in chemical 
potential, which, in turn, induces a great reverse pressure gradient, creating 
the environment of isothermal fluid distribution everywhere throughout the 
cylinder and homogenous distribution of superfluid density concentration. This 
illustration of the possible separation of superfluid from normal fluid of helium 
II means that there is in reality, no way for anyone to achieve the separation 
of the superfluid from the normal fluid of helium II. In other words, in 
considering the dynamical behavior of helium II in a large rotating cylinder, a 
mixture of the superfluid and the normal fluid without separation of the two 
fluids is accounted for in the model computation. Density concentration of 
superfluid is a function of temperature, which is also true for the surface 
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tension and viscous coefficient for helium II (Wilks, 1967; Hoare et al., 1961; 
Hung, 1990; Hung and Lee, 1992). In fact, the experiments carried out by Mason 
et al. (1978) showed that the classical fluid mechanics theory is applicable for 
cryogenic liquid helium in a large container. In this study, the theory of 
viscous Newtonian fluids is employed with modification of transport coefficients 

adjusted by normal and superfluid density concentration which is a function of 
temperature . 

II. Functions of Scientific Observation and Spacecraft Motions 
The GP-B spacecraft is a sun synchronous Earth satellite orbiting at 650 
km altitude directly over the poles. The functions of scientific observation for 
the GP-B spacecraft and its motions are illustrated as follows: 

The GP-B spacecraft is a relativity gyroscope experiment to test two 
extraordinary, universified predictions of Albert Einstein's general theory of 
relativity (Wilkinson et al., 1986; Stanford Relativity Gyroscope Experiment, 
1986). By using gyroscopes (ones with electrically supported spheres, spinning 
in a vacuum and others which utilize the spins of atomic nuclei, circulating 
sound waves, and even circulating laser beams), the GP-B measures two distinct 
space- time processes, frame (fraggjpg and the geodetic effect , which gradually 
changes its directions of spin. In these gyroscopes, the underlying principle 
is that rotating systems, free from disturbing forces, should remain pointing in 
the same direction in space. 

The physical meaning of "the same direction in space” can be explained from 
the differences in the Newtonian and Einsteinian concepts of the universe. For 
Newton, the answer was easy, in that space and time were absolutes. A perfect 
gyroscope set spinning and pointed at a star would stay aligned forever. Not so 
for Einstein, for whom space and time is warped. A gyroscope orbiting the Earth 
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finds two distinct space- time processes, frame dragging (rotation of space -time) 
and the geodeti c effect (curvature of space-time), which will be measured by the 
GP-B. 

As to frame -dragging for disturbances in the rotation of space-time, 
calculations based on a gyroscope in polar orbit at 650 km altitude should turn 
with the Earth through an angle amounting, after one year, to 42 milliarc-seconds 
(1 milliarc- second — 10 3 arc -second; 3.6 x 10 3 arc -seconds — 1°, angular degree) 
(Everitt, 1990) . As to the geodetic effect for disturbances in the curvature of 
space -time, the gyroscope's motion for rotation in the orbit plane at a 650 km 
altitude, from relativistic correction to the complicated motion of the Earth- 
Moon system around the Sun, should turn 6,600 milliarc-seconds per year (Everitt, 
1990) . 

f 

Furthermore, one of the GP-B's most important tasks is to investigate the 
gravitational action of moving matter in the universe. In other words, matter 
moving through space -time can be thought of as creating a new force - 
gravitomagnetism, which is similar to electromagnetic force created by an 
electric charge moving through magnetic fields (Misner et al., 1973; Everitt, 
1990). The frame -dragging measurement detects this force and fixes Its scale. 

In order to insure these extremely precise and accurate measurements and 
gyroscope operations, near zero temperature is required for mechanical stability 
of the instrument, preservation of the lead bag magnetic shield, shielding the 
gyroscopes against nongravitation disturbances, and for reading their direction 
of spin. Near zero temperature cryogenic liquid helium II (1.8°K) has been 
chosen to serve this purpose (Mason et al., 1978; Hung, 1990). 

The GP-B instrument comprises four gyroscopes and a reference telescope 
sighted on Rigel, a bright star in Orion. In polar orbit, with gyro spin 
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directions also pointing toward Rigel, the frame dragging and geodetic effects 
come out at right angles, each gyroscope measuring both. The GP-B attitude 
controller with a fine pointing system inside the dewar probe will constantly 
assure the operation of the spacecraft within the range of the scope of each 
experiment through the proper control of thrust. 

Thrust comes freely from the dewar' s boil-off helium vapor, in contrast to 
a traditional gas jet control system which uses on-off valves to save gas. Here, 
the copious supply of vapor allows the use of proportional thrusters : pairs of 
opposed nozzles adjusted continuously for much smoother control. Six thruster- 
pairs, suitably arrange, cover all three control functions with redundancy in 
case of failure. To take advantage of the characteristics of viscous -free 
superfluid helium II under X transition point, a porous plug is used to control 
the helium boil -off and to produce the fountain pressure to serve as their thrust 
jet for the spacecraft propulsion system. 

To comprise these functions of scientific observation, GP-B stores its 
gyroscopes, telescope, probe mass and others in the center core of the dewar 
probe surrounded by the cryogenic helium II liquid. The dewar container of the 
GP-B is spinning with a rotating rate of 0.1 rpm during normal operation. As the 
telescope is constantly sighted on Rigel, gyro spin directions also pointing 
toward Rigel during the spacecraft moving around the polar orbit. The GP-B with 
its rotating dewar, there is no slew motion involved in the spacecraft motion. 

III. Basic Charcteristics of Gravity Gradient 
and Gravity Jitter Accelerations 

Any fluid element inside the on-orbit spacecraft fluid system is subject 
to the acceleration that arises from the gravity gradient of the Earth 
(Avduyevsky, 1984; Forward, 1982; Misner et al., 1973; Hung and Pan, 1993; Hung 
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et al., 1993 a,b; 1994 a,b). Once the spacecraft orbit is fixed, the orbit 
period is determined and the basic structure of the gravity gradient acceleration 
also can be calculated. However, gravity gradient acceleration acting on each 
fluid element inside the on-orbit spacecraft fluid system is different dependent 
upon the distance of the location of the fluid element to the geometrical center 
of the spacecraft and its direction toward the location of the center of the 
Earth. This acceleration can only be calculated based on the non- inertia frame 
of spacecraft bound coordinate. Thus, the coordinate system shall be transformed 
from ordinary inertia frame coordinate to non- inertia coordinate. 

(A) Orbit Motion of Spacecraft 

Let us consider the case of the GP-B spacecraft, which is the Earth 
satellite orbiting at 650 km altitude directly over the poles, the orbit period, 
r a can be computed from following expression: 


, 3/2 


T„=2w- 


R e 9o 2 


where R E denotes the radius of Earth (- 6373 km); R c , the radius of the circular 
orbit (- R E + h - 7023 km); h, orbit altitude (- 650 km); and g,,, Earth gravity 
acceleration (- 9.81 m /s 2 ). For the case of the GP-B spacecraft, the orbit 
period r Q - 97.6 min, and orbit rate n - Zk /t 0 - 1.07 x 10" 3 rad /«- 

As the spacecraft is orbiting around the Earth, the azimuth angle of the 
Earth, toward the location of the spacecraft geometric center varies with 
respect to time. At time t — 0, the rolling axis of the spacecraft is aligned 
with the radial direction of the Earth's center to the spacecraft geometric 
center. Assuming the spacecraft rolling axis is linearly turning around 0° to 
360° in the orbit period, r Q , of the spacecraft when the spacecraft is orbiting 
around the Earth. This is particularly true for the case of the GP-B spacecraft. 
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Without the spacecraft slew motion, the azimuth angle (^ 0 ) can be defined as 



where r Q is the spacecraft orbit period [defined in Equation (3-1)]; and t is the 
time measured from the instant when the direction of the spacecraft rolling axis 
is aligned with the radial direction of Che spacecraft geometric center to the 
center of the Earth. 

(B) Spinning and Slew Motions of Spacecraft 
For the purpose to carry out wide-range observations, some scientific 
spacecraft requires slew motion with respect to the mass center of the 
spacecraft. For the case of the spacecraft slew motion, azimuth angle, shown in 
Equation (3-2), shall be modified through the coordinate transformation of slew 
motion when the spacecraft is orbiting around the Earth. 

Let us assume that the slew motion starts with the center located at the 
mass center of the spacecraft. Let us choose cartesian coordinate (x" , y" , z") 
with z"-axis along the axis of the dewar container (see Figure 1). At time t - 
0, the radial vector r c from the center of the spacecraft to the center of the 
Earth lies on the x"-z" plane of the cartesian coordinate chosen. The azimuth 
angle ^ is defined as the angle between the radial vector r fl and the z"-axis. 
Rotation matrices for spinning and/or slew motions along the x"-, y"- and z"-axes 
can be expressed as 

10 0 [coso) y t 0 -sinc*> y tl coscii x t sino),t 0 

0 cosc*) x t sinco x t 0 10 -sinu x t cosco^t 0 

0 -sin <i> x t cosu x £ sinca y t 0 cos<o y t o 0 1 

respectively. Here, « x , cjy and w, denote angular velocity of slew and/or 
spinning motions along the x"-, y"- and z"-axes, respectively. Radial vector r c 
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in cartesian coordinate without slew and spinning motion is 


f co = [siniJ; 2o/ 0, -costjrjJ (3-3) 

With an execution of spinning motion along the z"-axis only, radial vector r c 
becomes 


coso^t sin<i) r fc 0 





sin**. 

-sino) r c coso ) x c 0 


0 

0 0 1 


-cosi|r go 


= [3int So coso>,t, -sinilr^sinw.fc, -cosily] ( 3 - 4 ) 

With an execution of slew motion along the y"-axis only, radial vector r e becomes 


coso> y t 0 -sino) y t 


- 



sin**, 

0 l o 


0 

sinoj y t 0 coso>yt 


-cos^ 


= (sin(i|r 2o + a) y t) , 0, -cos (iJr^+cOyt) ] (3-5) 

With an operation of slew motion along the x"-axis only, radial vector r 0 becomes 


■ c-x 


1 

0 

0 


0 0 
cosu^t sinoj^t 
-sinoj^t coso > x t 


sin+ 2o 

0 

-cos**. 


= (sintlr*,, -costas inco x t, -cosilr^coso^t] (3-6) 

In other words, radial vector r e will be modified from the mathematical 
expression shown in Equation (3-3) to (3-4), (3-5) and (3-6) for the slew and/or 
spinning motions along the z"-, y"-, and x"-axes alone, respectively. In 
P art: f cu l ar > for the case of slew motion along the y"-axis, comparison between 
Equations (3-3) and (3-5), it shows that the azimuth angle will be modified as 
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so + w y C 


(3-7) 


For the successive operations of the spacecraft from spinning motion along 
the z"-axis, then slew motion along the y"-axis, and then slew motion along the 
x"-axis, radial vector r e results 


10 1 


cosWyfc 0 -sinu y t 

0 coso> x t sinu x fc 


0 10 

0 -sino ) x c cos o x c 


sincjyC 0 cosojyt 


cost > r t sino) r t 0 
• -sina> z C costa* 0 
0 0 1 

In addition to the modification of the azimuth angle made by the spacecraft 
slew motion through the formulation of coordinate transformation, shown in 
Equation (3-3) to (3-8), accelerations are also induced to activate on the fluid 
mass in the dewar container. Accelerations acting on the fluid particle in the 
dewar induced by the slew motion of the spacecraft with the coordinate fixed at 

the spacecraft center of the mass is as follows (see Figure 1) : 

»# 

R p * &x(&xR p ) + axif p + (3-9) 

where denotes the position vector of the fluid particle in the dewar container 
relative to the body frame of the spacecraft; c5*, angular velocity of the 
spacecraft body frame; a, angular acceleration of the spacecraft body frame; and 
v, velocity of the fluid particle relative to the spacecraft body frame. 

As we indicated earlier, let us assume that the slew motion starts with the 
center located at the spacecraft mass center, cartesian coordinate (x 1 9 , j* 9 , 
z* 9 ) is chosen with origin located at the spacecraft mass center. Let us also 
assume that x' ' -z ' 9 plane intersects the center of Earth and the spacecraft mass 


sinijr^ 

0 

-cosijr*, 


(3-8) 
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center. In other words, azimuth angle of Earth toward the spacecraft mass center 

lies in the x''-z'' plane. Slew motion is along both the x''- and y M - 

* « 

coordinates. Thus, w s - (w ax , u> ay , 0) and cT % - (a #x , a^, 0), due to slew motion 
becomes 


•r 

^p, slaw 



Ry 


+ «^r + 2w *y v '. 

3 

Ry H 

3 

-^(“sx^y-^sy) - a 3x R z~ 2(j3 sx V z 


R z n 

slaw 

+<a *y 2 > * +2(« 4Br v y -o Jr vr Jt ) 


(C) Coupling for the Accelerations of Spinning and Slew Motion of Spacecraft 
For some particular reasons required in the scientific spacecraft, it might 
be faced with the situation that both spinning and slew motions are needed 
simultaneously. To encounter this case, the following formulations are made to 
deal with coupling for the accelerations of spinning and slew motion of the 
spacecraft: 


cT 

p, and spinning 



slaw and spinning 


+ «sy*r + 2 « 5 y V * 
MxXy-UsyRx) - a sx R z- 2 ^sx V z 
-*z< w x 2+w «y a > + (<t ax R y -a ay R x ) + 2(u> ax v y -u sy v x ) 


J slaw 


- ( U **x - ««**) <*>r ‘ <M y - 2 <*z V yl 


(3-11) 


spinning and coupling 


where w x and w z denote angular velocity and angular acceleration, respectively, 
of spacecraft spinning motion along the z-axis. 

For the case of the GP-B spacecraft, there is no slew motion and the 
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spinning is the only acceleration acting on the spacecraft fluid system. 
Acceleration due to spacecraft spinning motion becomes 


*p, spinning 


R, 


spinning 


-R x u> z 2 -R y C} s -2u> z v y 
-J? y u r 2 +J? x « z + 2u z v x 


(3-12) 


spinning 


To convert the expression of Equation (3-12) in cartesian coordinate to 
cylindrical coordinate, by using the relationships of (R*, Ry) - (rcostf, rsin*) 
and (v x , v y ) - (u* cos* - u„ sin*, u* sin* + u* cos*). Equation (3-12) becomes 





X 


-rcos 0 o) z 2 -rsin 06 r -2 (u r sin 0 + UgCosO) w x 

spinning * 

if y 

n 


-rsin0o r 2 +rcos0<i x + 2 (u r cos0 -i^sinO) <i> x 



spinning 

0 


(3-13) 

spinning 


and 


*p.. 





spinning ~ 

r? 

s 



spinning 


R x cosQ + RySinB 
-R x sinQ + RyCosB 


-r« x 2 -2ue« x 

rt!) x + 2u r oi x 

0 


z J spinning 

Accelerations induced by spacecraft spinning motion alone becomes 


(3-14) 

spinning 




kl 

m 


ro) x 2 + 2 uj(i) x 

a a 

= — 


= 

-ru x -2u r o> x 


spinning 


spinning 

0 


(3-15) 


spinning 


(D) Gravity Gradient Acceleration 

The gravity gradient acceleration acting on the fluid mass of spacecraft 
can be shown as 


£„=n 2 [3 (f c -<3) t c -d] (3-16) 

where a gg denotes gravity gradient acceleration vector; <1, the vector (not a unit 
vector) from the fluid element to the spacecraft geometric center; £ e , a unit 
vector from the spacecraft geometric center to the center of the Earth; and n. 
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the orbit rate (see Figure 1). 

It is assumed that the gravity gradient exerted on the geometrical center 
of the spacecraft orbiting around the Earth on its specified orbit is zero. In 
other words, all the gravity acceleration exerted on the spacecraft is nothing 
but the gravity gradient acceleration which is defined in Equation (3-16). In 
this study, we are interested in investigating how gravity gradient acceleration 
affects the dynamical behaviors of cryogenic fluid elements of helium. 

For the convenience of mathematical calculation, let us describe all the 
parameters involved in Equation (3-16) in terms of cartesian coordinates. In 
order to match with the computer simulation, mathematical derivation are 
considered in the first quadrant. Figure 1 illustrates the geometrical 
relationship of the parameters shown in Equation (3-16). 

Let us consider the fluid element of interests, m, located at (r, 9 , z) in 
cylindrical coordinates and at (x, y, z) in cartesian coordinates. The origin 
of the two coordinate systems is located at the center bottom of the dewar tank. 
The slew and/or spinning motions, mentioned earlier, are executed at the 
spacecraft mass center with cartesian coordinate (x", y", z"). The geometry 
center or the spacecraft mass center is located at z - 4 . As |d| (not an unit 
vector) is much smaller than the distance between the location of the GP-B 
spacecraft geometric center to the center of the Earth, r e (an unit vector) 
through the GP-B geometric center and r c (an unit vector) through the fluid 
element, m, is basically the same. Assume that vector r e lies in the x-z plane 
of the cartesian coordinate. 

Radial vector r c with the modification of slew and/or spinning motions 
along the x"-, y"-, z"-axes have been derived in Equations (3-3) to (3-8). Based 
on the relationship between coordinates (x, y, z) and (x", y", z") 
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(3-17) 


X 


1 0 0 


x" 

y 

= 

0 10 


y" 

z-L^ 




_// 

c 


o 

O 


z 


Vector 8 in (x, y, z) coordinate becomes 


d ~ [-rcos0, -rsin0, -(z-L c ) ] (3-18) 

Substituting Equations (3-4) and (3-18) in (3-16), non-inertia frame 

expression of gravity gradient acceleration with spinning motion in z-axis 
becomes 




• n 2 


3 (~x*ia^ fc cox {0 +* g t) ♦ (x-L^ cox^J lin^^coxca^t + rcoxS 
-3 C-rBin*j*co«(e«.tt,t) * «intj*xin«,C + r*in6 

C~*’*in^ JD cax (6 + <*,£) ♦ (x-£ c ) cox^jJ coxf (x-L^J 


(3-19) 


Substituting Equations (3-5) and (3-18) in (3-16), non- inertia frame 


expressions of gravity gradient acceleration with slew motion in y-axis becomes 


a 99.x 


3 [-rsinijfcosO +cos^r ( z-L c ) ] sin* +rcos0 

a ss,y 

= n 2 

rsin0 

a 97- *. 

slaw in y-axis 

-3 [-rsin*cos0 +cos* (z-L c ) ] cos* + ( z~L c ) 


where - ^ + Wyt. 

Substituting Equations (3-6) and (3-18) in (3-16), non- inertia frame 
expressions of gravity gradient acceleration with slew motion in x-axis becomes 


fa J 


m.x 


*-■* 

« a 2 




3l-rcoa6(iaf ll *eof f B [riinu ( ttint*cosu ( t (z-U ])iinf A *reoi( 
Oi-rcoxdxinf^+coxtjJ 1 »in«,t»ind*cox« x t( ] teoat^sintt-t+rxinO 
-3(-rcos0xinf j*+cox [rxin<i> x eaind+coxw x e ( x-L*) J kjoxf ^cox^e* ( x-Z*) 


(3-21) 


The relationship for the coordinate transformation from cartesian to 


cylindrical coordinates for any vector F (such as velocity or force vectors ) in 
non- inertia frame of spacecraft bound coordinate can be shown as 
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(3-22) 


\f l [ cos0 sin0 ol f 1 

x r x 

F e = -sin0 cos0 0 F y 

F z\ 0 0 1J [ F * 

Thus, the gravity gradient acceleration located at (r, 6 , z) can be computed from 
that located at (x, y, z) , shown in Equations (3-18) to (3-21), from the 
following relation: 


'^wl f 0080 sin0 °1 pW 

£gg = *99.0 = -Sin0 COS0 0 &gg ,y 

*99. z 0 0 1 * 99 , Z. 


(3-23) 


(E) Jitter Accelerations 

In addition to gravity gradient acceleration acting on the fluid element 
of on-orbit spacecraft fluid systems, there is another acceleration of gravity 
jitter also exerted forces on the fluid systems. The sources of residual 
acceleration of gravity jitter range from slew motion of spacecraft, atmospheric 
drag on the spacecraft, background gravity, spacecraft attitude motions arising 
from machinery vibrations , thruster firings, crew motion, etc., are also capable 
for the excitation of slosh waves in spacecraft fluid systems (Kamotani et al . , 
1981; Hung and Shyu, 1991 a,b,c; 1992 a,b,c; 1993 a,b; 1994 Hung et al., 1992 
a ,b ; 1994). 

Among all of the varieties of jitter accelerations listed, accelerations 
induced by slew motion of the spacecraft dominate the forces activated on the 
spacecraft fluid systems. Two coordinate systems (cylindrical and cartesian) 
chosen in this study are (r, 0, z) with corresponding velocity components (u r , 
u fl , u 2 ) for cylindrical, and (x, y, z) with corresponding velocity components 
(u*, Uy, u z ) for cartesian coordinates. The origin of these two coordinates are 
located at the central bottom of the dewar tank, as shown in Figure 1. The 
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spacecraft center of mass, or the geometric center of the spacecraft is located 
at (x c , y c , z e ) - (0, 0, L c ) . The relationships of the coordinate, velocity and 
the force between cartesian and cylindrical coordinates are 


x 


COS© 0 0 


z 

y 

= 

sin0 0 0 


0 

z 


0 0 1. 


z 


N 


cos0 -sin© 0 


'«r 

“r 

= 

sin0 cos© 0 


«• 



0 0 1 




\ F r] 


cos0 sin0 O' 


F 

A X 

F * 

= 

-sin0 cos0 0 


F y 

F z\ 


0 0 1 




(3-26) 


In the derivation of acceleration induced by the slew motion of spacecraft, 
the coordinate system (x" , y" , z") is fixed at the spacecraft center of the mass. 
The relationships of the coordinate, velocity and acceleration between 
expressions with the origin located at the spacecraft center of the mass (x" , y" , 
z w ) and origin located at the center bottom of the dewar tank (x, y, z) are 
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COS0 0 0 


r 

y 

s 

sin0 0 0 


0 

z-L c \ 


0 0 1. 


. Z ~ L c. 


(3-27) 


[v x 1 
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-1 

l*rl 

i? y 
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F y 
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-sin0 cos© 0 
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SitfV 

0 0 1 


F z 


(3-29) 
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(3-30) 


pv 


' cos 0 sin0 0 


= 

-sin0 cos© 0 

F z 

slow 

0 0 1. 



slow 


A detailed expression of [R*, Ry, R,],^ are shown in Equation (3-10) of 
this report. Jitter acceleration is a summation of acceleration induced by slew 
motion and others, such as atmospheric drag on the spacecraft, spacecraft 
attitude motions arising from machinery vibration, thruster firing, crew motion, 
etc. Thus, jitter acceleration can be expressed as 


- 

'«W 

+ 

a <ri,e 

s — 

pv 

F* 


\F e ] 

F 9 

{ l + -i sin (2*^t) } 



alow 

a n.z. 

' others 

Fg\ 

slow 

w 

2 

others 



‘cos0 sin0 O' 


r f i 

X 


\ F r] 


= - 

-sin0 cos0 0 


Fy 

- 

F 9 

{ 1 + -^-sin(2it.ffc) } (3-31) 


0 0 1 


w 

slow 

w 

2 

others 


where f is the jitter frequency (Hz) imposed on the fluid systems of the 
spacecraft . 

IV. Non- Inertia Frame Mathematical Formulation of Fundamental Equations 
Dynamical behavior of fluid elements inside the on-orbit spacecraft fluid 
systems are strongly modified by the gravity gradient and gravity jitter 
accelerations . In order to accommodate the impact of gravity gradient 
acceleration, in particular, on the on-orbit fluid motion, one has to consider 
non- inertia frame of the spacecraft bound coordinate rather than adopting inertia 
frame coordinate used in ordinary fluid mechanics formulation. 

Consider a closed circular cylindrical dewar of radius, a, with height, L, 
which is partially filled with cryogenic liquid helium, and rest of the ullage 
i- s filled with a helium vapor. Angular velocity of rotating cylinder is w. 
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Density and viscosity of liquid helium and helium vapor are p L , /i L , p v , and 
respectively (Mason, 1978). Let us use cylindrical coordinates (r, 9 t z) , with 
corresponding velocity components (u, v, w) , and corresponding residual gravity 
acceleration, such as gravity gradient components (a gg r , a gg a gg 2 ) anc * gravity 
jitter components (a g j r , a g j^, a gj z ) . In the derivation of the governing 
equations, accelerations induced by the spinning motion of the spacecraft is 
included in the formulation. The rest of the acceleration such as slew motion, 
atmospheric drag on the spacecraft, spacecraft attitude motions arising from 
machinery vibrations, thruster firing and others, are included in the jitter 
acceleration, shown in Equation (3-31). The governing equations for non-inertia 
frame of spacecraft bound coordinates can be shown as follows: 

(A) Continuity Equation 


4 4 (r U ) + 4!£*J^0 


(4-1) 


(B) Momentum Equations 


_ , du . „ du . v du v 2 

p ( ac +u a7*71e-7 


) =- +2p<ii r v+p (a. 


gj,r J> ’ a gg, 


,,) + P rW z 2 


+p (^u- — 
i 2 


2 dv v 

r 2 ae ; 


(4-2) 


P 


, dv u dv _v dv _uv 
at dz r 30 r 


dv > 

+ V- 3 -) 

dz 


1 dp 

z ae 


-2po*u+p (a^ /0 + a ?ff/ a) ~P^* 


+ |i (^v-— — 


2 du 

r 2 ae 


) 


(4-3) 




where 
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(4-5) 


a )+ _LJL + . a 2 


r dr dr r 2 d6 2 dz 2 
In these formulations, 2w z v and 2o> z u are the Coriolis acceleration, rw x 2 is 
the centrifugal acceleration, and rw z is the angular acceleration induced by the 
spinning motion of the spacecraft. 

In the computation of fluid forces, moment, viscous stress and angular 
momentum acting on the container wall of the spacecraft, one has to consider 
those forces and moment in the inertia frame rather than the non- inertia frame. 
In particular for the case of the spinning motion in the z-axis. To show an 
example, one has to transform those vectors from the non- Inertia frame to the 
inertia frame for the case of spinning motion in the z-axis. 


F r 

r X 


cost o g t -sino> z t 0 


Pel 

K 

3 

sinu z t costo z fc 0 


F y 

w 


0 0 1. 


w 


(4-6) 


where prime symbol denotes vectors in the inertia frame while those parameters 
without the prime symbol indicate vectors in the non- inertia frame. 

V. Initial and Boundary Conditions of 
Spacecraft Fluid System in Microgravity Environment 
Governing equations of the fluid motion in on-orbit spacecraft fluid 
systems for non- inertia frame of spacecraft bound coordinates have been 
illustrated in Equations (4-1) to (4-6) . These equations shall be combined with 
the characteris tics of gravity gradient and gravity jitter accelerations as 
formulated in Equations (3-1) to (3-31). Initial and boundary conditions shall 
be introduced to accommodate solving fluid motion in on-orbit spacecraft fluid 
system for non-inertia frame coordinate (Hung et al. f 1990 a,b,c; 1991 


a,b,c,d,e,f,g,h,i,j). 
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Lee the profile of the interface between gaseous and liquid fluids be given 

by: 

(5-1) 

The initial condition of the profile of the interface between gaseous and 
liquid fluids at t - t„ is assigned explicitly, and is given by: 

'K t “t 0 , r, 9 , z) -0 (5-2) 

A set of boundary conditions has to be supplied for solving the equations. 
These initial interface profiles used in this study have been given explicitly 
through the steady state computations made by Hung and Leslie (1988) and Hung et 
al (1989 a,b,c,d) which were checked by the experiments carried out by Leslie 
(1985). These boundary conditions are as follows: 

(1) Along the container wall, the following three boundary conditions 
apply: 

(a) Interface between solid and liquid fluid: No -penetration and no-slip 

conditions assure that both the tangential and the normal components 
of the liquid velocity along the solid walls will vanish. 

(b) Interface between solid and gaseous (vapor) fluid: Similar no- 

penetration and no- slip conditions as that shown for interface between 
solid and liquid fluid will apply. 

(c) At the location of solid- liquid-gaseous (vapor) three phases 
interfscs . No-penetration, but not no-slip condition apply. This 
will assure that normal components of liquid and vapor velocities 
along the solid wall vanish, and allow a slipping flow of liquid and 
vapor fluids along the solid wall at three phase interface location. 
The velocity of slipping flow at this location is governed by the 
adhesive forces between fluids (liquid and gaseous) and solid walls. 
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Also, at this location of three phase interface, a constant contact 
angle is present in which the behaviors of wet or dry contacts are 
determined by Coulomb interaction between the fluids (liquid and 
vapor) and the surface phenomena (material and roughness) of solid 
walls . 

(2) Along the interface between the liquid and gaseous 

fluids, the following two conditions apply: 

(a) Kinematic surface boundary condition: The liquid (or gaseous) surface 
moves with the liquid (or gas) which implies 

^5r=°,or 

Dt 

(s-3) 

at or r dO oz 
on r\ ( t=t a , r,0, z) 

(b) Interface stress condition: Across the liquid-vapor interface, the 

stress must be continuous. Based on Landu and Lifshitz (1959) f the 
stress across the liquid-vapor interface can be expressed as 

lP 9 -P L )n r Ux a ) 0 -(T„) L -\n r oi-i- * -ji-ln, (5-4) 

where R x and R 2 are the radius of curvatures of two major axes at the point of 
interests on the surface of the liquid- vapor interface. 

The expressions of radius of curvatures R x and R 2 in cylindrical 
coordinates from differential geometry can be shown as 

1 1 __ EN-2FM+GL /k_r\ 

(5 5) 

i?i *2 EG-F 2 

where the relationship of cartesian and cylindrical coordinates for the curved 
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surface of liquid-vapor interface is 
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rcos0 

y 

3 

rsin0 

z 


H(C',r,B) 


Here, the configuration of the liquid-vapor interface is z - H (t - 


(5-6) 


t., r, 9) 


e - (t r * (t r * (ii)' 


(5-7) 


F = dx + + dz 8g 

8r 80 8r 80 dr dQ 


(5-8) 


where , 
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{EG-F 2 ) 2 
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(EG-F 2 ) 2 
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Z& 
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(5-9) 


(5-10) 


(5-11) 




(5-12) 


t*®/ yL Z& = -^[*,y,z] 


(5-13) 
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and, 




a 2 

// ..!> // 



■^rr Xrr ^rr 


dr 2 

** y* 2 % 

s 

a 2 

drae 

y* 


a 2 

ae 2 . 


Simplifying Equations (5-7) to (5-11), one obtains 

EG-F 2 = r 2 (l + H 2 + —a 2 ) - Z 2 D 2 

r 2 

E - 1 + H 2 , F - H G = r 2 * Hf 


(5-14) 


D = 


(1 ♦ H 2 * JjV) 

r 2 


i 

7 


Ln Os.. Mm ilB. * Z JL rt) „ _ (X 2 H r + 

0 ru ' ^ rD ' 


and 


Substituting these relations to Equation (5-4), the radius of the curvature 
on the curved surface of the configuration of liquid-vapor interface in 
cylindrical coordinates can be expressed as follows! 


_1_ 




dr D 


ae in ' 1 


Here, in Equations (5-4) and (5-15) 






2 tot, 

3 ax* 




(5-15) 


is the viscous stress tensor; ji, the viscous coefficient of the first kind; f. 
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the viscous coefficient of the second kind; P, the pressure; a, the surface 
tension of the liquid-vapor interface; and n i , the unit vector normal to the 
interface; and 6^, the Kronecker's delta function. Also, subscripts G and L 
denote conditions at gaseous and liquids fluids, respectively, across the liquid- 
vapor interface. 

The fluid stresses across the liquid-vapor interface can be decomposed to 
the components normal (n A , a unit vector) and tangential (t i( a unit vector) to 
the interface. For the component tangential to the interface, one can take a dot 
product of a unit vector tangential to the interface, t it to Equations (5-4) and 
(5-15), which leads to 

t ( ti n j) ] t ( t ij t-ifij) 1 q (5 — 16) 

since n^ - 0. 

For the component normal to the interface, one can also take a dot product 
of a unit vector normal to the interface, n i( to Equations (5-4) and (5-15), 
which leads to 


Po ~ 


P L - I <w,».- ♦ 4<-£n 


* 


For components normal to the interface along the (r, 6, z) directions in 


cylindrical coordinates can be obtained by taking dot products of n r> n*. n* 
separately to Equations (5-4) and (5-15), which are expressed as 



' n z 


( X ) G - ( X rj^j ) L 

(Pa ~ Pi) 

n z 

L 2 J 


G “ ( X Gj n j) L 
( T Zj^j ^ G “ ^ X Zj**j ) l 


where (n r , n d , n z ) is the unit vector normal to the interface in cylindrical 
coordinates (r, 6, z). 
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(5-18) 


= __o [ d { rHj L) + 3 (A)j 
r dr D 30 rD 




For a special case of axial symmetry, component normal to the interface, 
shown in Equation (5-17), can be simplified and becomes 


P a ~ P L ~ ( x lj n i n l) <3 * ^ T i J n l n j) l 


^ -] (5-19) 

r dr 1 

(i + 4> 2 ) 2 

since d/dd - 0, H r - dH/dr - dz/dr - H$ - 0 and D - (1 + <fr 2 ) in for the case 
of axial symmetry. 

VI. Characteristics of Slosh Wave Induced Fluctuations 
in Fluid Moment and Angular Momentum 

Slosh wave Induced fluctuations in the fluid system of the rotating dewar 
introduce time -dependent disturbances In moment and angular momentum of 
spacecraft fluid system. In this study, there are induced angular velocities 
along the yawing, pitching and rolling axes due to the fluid motion inside 
rotating container. These angular velocities in yawing, pitching and rolling 
axes, caused by the fluid flows in a partially liquid-filled container, readjust 
the angular velocity in rolling axis. 

In order to accommodate the spacecraft dynamics of yawing, pitching and 
rolling, cylindrical coordinates (shown in Figure 1) of rotating container is 
transformed into cartesian coordinates based on (x, y, 2 ) - (rco sB , rslni*, z) 
with corresponding velocity components (V x , V y , V x ) - (ucos* - vsin* , us In* + 
vcos*, w) . If spacecraft is rotated with respect to mass center at (r c , 9 ct z c ) 
in cylindrical coordinates , location of mass center in cartesian coordinates 
becomes (x c , y c , z c ) - (r c cos* c , r c sin0 c , z c ). Induced angular velocities (w*, wy. 
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tj x ) in cartesian coordinates becomes 


where 


fc3 1 


*** -^xy **z 


\ V *\ 

a 

y 

= 

Kyx Kyy ~^yz 


V y 



~^zx ^zy K zx 


[Vz\ 





{Z-Z a ) 


t (y-y c ) 2 +(z-z c ) 2 ] -1 
[ (x-x c ) 2 +(z-z c ) 2 ] -1 


(6-1) 
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As the velocity components are given by 


v x 


cos© -sin0 0 

V y 


sin0 cos0 0 

rJ 


0 0 1. 


U 

V , 

V 


(6-2) 


the relationship between the components of induced angular velocity and flow 
velocity in cylindrical coordinates can be expressed in the following 
formulation: 


where 
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For the case of the GP-B, the axis of rotation is always fixed at the point 
of proof mass which is located at the geometric center of the dewar at (x,, y c ,z c 
— (0, 0, Le) where 4 . - 1/2 L and L is the height of the dewar (see Figure 1) . 
By using the computed results of induced angular velocity shown in Equation (6- 
3), one can compute the angular momentum (H*, Hy, H,) as follows: 
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and 
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The moment of spacecraft can be computed from the time rate of change of 
the angular momentum, i.e. f 
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(6-5) 


where - (w x , Wy, w x ) denotes pitching, yawing and rolling angular velocities 
of spacecraft in inertia frame. 

VII. Mathematical Formulation of Fluid Stresses 
and Moment Fluctuations Due to Slosh Waves 
For the purpose of considering large amplitude slosh wave activated fluid 
stresses exerted on the solid walls of the dewar, the fluid stresses are 
decomposed into the tangential and normal components to the walls which can be 
expresses as follows : 



(7-1) 


n =p5 




(7-2) 


where n t denotes the tangential component of fluid stresses; Qq, the normal 
component of fluid stresses; P, the thermodynamic pressure; u^, fluid velocity 
in a direction; t a , unit vector tangential to the wall; n^, unit vector normal 
to the wall; /i, the molecular viscosity coefficient of fluid; and S^, the Dirac 
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delta function. Subscripts a and 0 imply the directions of flow fields. 


Figures 2(A) and 2(B) show the geometry of GP-B dewar propellant tank in 
both r-z and r-0 planes, respectively. In order to make the computation of fluid 
stresses match the geometry of the dewar tank, mathematical formulations have 
been divided into three sections: (A) Top wall (dome) section, (B) Bottom wall 

(dome) section, (C) Cylindrical section (including probe section of the inner 
wall of dewar) and (D) Baffle plates section. There are several plates of baffle 
inserted in the dewar [see Figure 3(A) and 3(B) for baffle board illustration] . 
Baffles with the shape of hollow circular plate with an inner radius R lt and an 
outer radius Rj are installed along the probe column of the dewar located at z 

- 1^, where i - 1, 2 n plates and the thickness of each plate is d. Figure 

4 shows the GP-B dewar container equipped with probe and baffle boards. 

(A) Top Wall (Dome) Section: 



(n e >£! +!f)cos24> 

(7-3) 

®-> g If ^cos^lf^sirKt.] 

(7-4) 



(7-5) 

(B) Bottom Wall (Dome) 

Section: 
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(7-8) 
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(C) Cylindrical Section (including probe section of the inner wall of dewar) : 


(D) Baffle Plates Section: 


<n c ) 'cylindrical**^ f + Jj) 

(7-9) 

(D t ) ^ ij3dricaJ =n|-i 
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Cylindrical* ^ 

(7-11) 

<nj r - r »n(iH + i^) 

c Baffle dz dr 

(7-12) 

/TT \ r-9 _ , 1 du dv. 

^Baffle - * a? ) 

(7-13) 

Baffle 3 

(7-14) 


where <0 is the azimuth angle of the dome; and (lit)*** and (II t ) r ~* denote tangential 
stresses in r-z and r-0 planes of the dewar, respectively. Velocity components 
in cylindrical coordinates of (r, 0, z) are shown as (u, v, w) . 

The stress distribution shown in Eqs. (7-3) to (7-14) can be integrated 
with respect to area and obtain the tangential and normal forces on top wall, 
bottom wall, and cylindrical sections of the dewar. 

(A) Top Wall (Dome) Sections: 
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(B) Bottom Wall (Dome) Section: 


'(Ft) 1 -* 

•// 

' (n t ) 

(F c ) r -* 

(n c )'-« 

(F a ) 

Bottom Well 

<n a ) ] 


J?3cos<|>d4>d0 


(7-16) 


\Bottam Wall 


(C) Cylindrical Section (including probe section of the inner wall of dewar) : 
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Outer Cylinder 


(D) Baffle Plate Section: 
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where and <f> denote the radius and the azimuth angle of the dome section, 
respectively; and and R^ are the radii of the inner and outer walls of the 
cylindrical section of the dewar (see Figs. 2 and 3). 

In order to accommodate the spacecraft dynamics of pitching, yawing and 
rolling, cylindrical coordinates (shown in Figs. 2 and 3) of the rotating 
container is transformed into cartesian coordinates based on (x, y, z) - (rcostf, 
rsintf , z) with corresponding velocity components (v x , v y , v 8 ) - (ucostf - vsintf, 
usin0 + vcostf , w) . For the case of the GP-B Spacecraft, the axis of rotation is 
always fixed at the point of proof mass which is located at the geometrical 
center of the dewar at (x c , y c , z e ) - (0. 0, L/2) , where L is the height of the 
dewar (see Figs. 2 and 3). To fulfill this goal, stress distributions , shown in 
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Eqs . (7-15) to (7-18), have to be recalculated in the (x, y, z) directions (see 
Figs. 2 and 3). 

M [(/-<"■> Top Wall * //< n ‘> Soceaa Well] i^COS<|>-sin<|>*COS0d<j>d0} 

Top Wall* //<“■> ««ji] K|cos<Min0d<t>d8} 

+{ [// <H a ) Top mL 1+// <&„) Socto. Mil] i?dcos 2 4>«cos0d<t>de} 

+ | £ f f ~ (2,) xn/!«rcyiind»r® + //<“»> lnn#rcyliiid#r cos ®] 

*£{ Iff < n => ’£n . cos6rdrda *[[ <n => £m. 3in0r * dB L,-2 • 

*5 | [//(n t )- fJ# cos8rdrdB ♦ (H e ) ^ J (7-19) 

F y ={ [[J- <n c ) ££ *.ji + / / (D c ) £«<» w.ii] i?dcos<t>*sin<|>*sin0d(t»de} 

+ < [//- (n c ) oue«r Cylinder 3in0 ♦ // (DJ cyi^COS©] * c *d0dZ ) 

Top Wa 11* //«■> 5£tt» Rdcos<|>*cos0d4>de} 

+{ Iff (D a ) rop hu* f f (D n ) Boccoa * 22 ] ^cos 2 4>-sin0d<t>d0} 

+{ Iff (n c )^. rcyiind#r cos0 + // (n n ) ^cyi^-rSine] * cl d0dz} 
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+ { [ff ^n) Outer Cylinder 003 ® + // Outer Cyl Inder 3 in® ]•£« ) 

| [//< n = )r s .« J . 3in9rcird8 * <n e > ««i. c:oserdrd9 ],. tl .| } 

*?{[[/ (n ' } r */™. 3in0xdrd0 * // (n ' ) r ^ fi . cos6xdrc * L tl .|} ( 7 ~ 2 °) 

F,=| j J j* (H c ) Top Wall*f J (E c ) soctoo Waii] i?dCOS 2 <|>d(|>d0| + U f J (E n ) Top wall 

(B. x ) lotnm mill Rjcoe4-ain4d><#)«|Jj’(II £ )ii c< yii^ r J? CJ ^d* ♦ j'J(n,)oBMrorji^r*e» < * d ^ ^ 

The moment of stress force acting on the dewar wall of the container can 
be computed from the cross product of moment arm, [which is the perpendicular 
distance from the location of proof mass of the spacecraft to the total forces] , 
and the total forces. Components of the moment can be computed from the 
following formulations : 

M x = UJ (n e ) Top wall* ff (n c ) toltom wall] i?dCos 3 <t>-sin0<*|>d0 

+[/ / (n a ) rop goccoa *, J2 ] ^cos 2 ^sin4>-sin0d«Dd0 

*<[//-« n ) out^r Cylind0z Sin0 + //«”■> Outer cyiiadar^®®®] 

* cniimtSiiatne&dzofl (H t ) oucax cyu^torRijainfliBdi] * ( (//" •“■**// »u] 

• (L c -L cf +i? d sin4))i?|cos<l)*sin^*sin0d^d0+ i//(n e) Top wall 
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*//«•> attorn mu! (L e -L d +i? (i sin<l>)i?^cos4>-cos0cf<l>d0 + t//<n 0 ) Top Wall 
* // (D a ) soccoffl WJ ] (L^-Lj+^inilr) J? d 2 cos 2 i|rsin0di|rd0 

{ [ J J (D c ) Ywarcyllad«r cos ^ + f J* (II n ) rnnar cylinder 3 ] ( z ~ L c) R o jd0dz 

* [ // «v Oue#rcyiind»rCOS0 + ^ ^ Out0i Cylinder ^ ^ L c ) i? c ^d0dz} 

| [-/ / <n «> t »r,i. 3in<>rdzcie *ff <n '> «m- C039r<±rcB ],. tl -| ' - L ^ L i' 4 ) } 

| [■// (n ' } /i . ainerdrd0 * // (IIt) w«. c ° 80 rdrc *L il .-f ( ~ l ** l ‘* 4 ] } ( 7_22 ) 
Af y =jJ / (H c ) toJ2 +/ / <n c ) Volt™ wnij](-C'c“ i 'd +s d s in4>)^dCos<|)*sin<J>-cos0d(|>d0 

+[//- (D c ) &J (D e ) £ e t to» Mil] (L c -L d +*dSin<|>) * J? d cos4>*sin0cS>dB 

+[[/ (E n ) T opt*u+ff (D.) Bottom mu] (Lc-^d3in<l>) ^d-cos 2 <j>-cos0dt>d0 
+ { UJ- (n c ) imer Cylinder 9 i- n Q+ff (E n ) lonar cylinder * COS ®] ' ( Z “^c)' R cl^® < ^ Z 
+ [// * (n c ) Oue# r cyiiad«r s in0 * // (n»> c»t«cyiW=° 3 «] (z-L c ) ^d0dz} 

-{ t J J (II C ) Top Hal 1*J f <n c ) Bottom Iftiil •f?dCOS 2 <J)*COS0d4>d0+ [J j (II n ) jbp 
+// (D a ) Boccomwaii] *d^os 2 <l>* 3 in<l>-cos 0 d<l>d 0 

+ J J (H c ) r^«rCylind«r^clCOs4)d4>d0 + J j (1I C ) Ouc« r cyiind*r^ci cos< t > dl|)d0} 
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{ [[/<“«> ^n. eo - #rdr< * * //<nc)^.ain0rdrd0] i _ t ^(-L c >L 1 -|)j 

{ (-// < n =) w/ I . CO80jrdr< * * // (n t J w«. 3inerdr< *L < .- ( ' t -' £ -‘ > T 1 } (7-23) 

M g ^J f - (n c ) ££ wall *f / (n t ) £Uam Waii ]i?dCOs 2 <|>-sin«|>-sine-cos0d(|>d0 
*jj'f (n e ) Top Bottom «iiiJi?d cos2( l > ‘ cos2 ® < ^ < J ><: ^® 

+(/ / (n D ) rop vai i + //(H n ) WJJ ]^cos 3 <|>-sin0*cos0d(|>d0 

+[// (n c ) ^ercyllndarc OS0+JJ* (D a ) Ina^r cyl index ^ ]^cl c ° 3 ® 

<*[//<“=> Outer cylinder^OsS * //<“•> OucarCyii n d.rSin0]i?J^OS0d0dz 

-{[//- (n c ) ££ Wali +J f (n c ) S^ con , * aJJ ] *3cos 2 <|>7jin4>^os0-sin0d(t>^ 

+[//- (n c ) £j wajj-j " f (n c ) £? ecoa Wa ii]i?3cos 2 <t>-sin 2 0<*t>d0 
+[[ f (H n ) rop » n + (n n ) floctaa MJJ ]^cos 3 «|)-cos0*sin0d(t»d0 

*[//-«=» CyUndar 3 i- n Q + f f (E n ) rnn«rCyJi/id«r COS ®]^-I 3 

*[//-«•» Ouear cyi iadar 3 i H0 + JJ (E n ) ouearcyiindar 003 © ] 7? ci Sin0d0dz} (7-24) 

where denotes the height of the dome (see Fig. 2). 

After the integration of Eqs. (7-19) thru (7-24), components of the moment 
can be rewritten as the following expression: 


\M X ] 


fV,-(L,-i c )F,l 

My 

3 

U,-L C )F X -LS, 

kJ 




(7-25) 


36 



where L*, Ly, and L, denote the components of the moment arm along x, y, and z 
axes , respectively . 

By using the relation of F a L* - 0, or F x L* + F y Ly + F* (L* - 1^) - 0, 
moment arms of the moment of fluid stress moment induced by the slosh wave 
excitation can be computed from the following relations: 

[ L X 1 . \FyM z -FM 

L y = i F Z M X ~F X M Z (7-26) 

L t -L c \ F*+F}+F* [FJiy-FyMx 

VIII. Methods of Numerical Simulation 

Detailed descriptions of the computational algorithm applicable to 
cryogenic fluid management under microgravity are also given in our earlier 
studies (Hung et al. # 1990 a,b,c). In this report, a full-scale GP-B and AXAF 
spacecraft propellant dewar tanks with a radius of 68 cm and a height of 145 cm 
will be used in the numerical simulation. The propellant tank is 80% filled with 
cryogenic liquid helium for the GP-B dewar while the rest of the ullage is filled 
with helium vapor. The temperature of cryogenic helium is 1.8 K. In this study 
the following data were used: liquid helium density - 0.146 g/cm 3 , helium vapor 

density - 0.00147 g/cm 3 , fluid pressure - 1.66 x 10 3 dyne/cm 2 , surface tension 
coefficient at the interface between liquid helium and helium vapor - 0.353 
dyne/cm, liquid helium viscosity coefficient - 9.61 x 10“ 3 cm 2 /s; and contact 
angle - 5°. The initial profiles of the liquid- vapor interface for the rotating 
dewar are determined from computations based on algorithms developed for the 
steady state formulation of microgravity fluid management (Hung et al., 1990 
a,b,c) . 

A staggered grid for the velocity components is used in this computer 
program. The method was developed by Harlow and Welch (1965) for their MAC 
(marker-and-cell) method of studying fluid flows along a free surface. The 
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finite difference method employed in this numerical study was the "Hybrid Scheme" 
developed by Spalding (1972). The formulation for this method is valid for any 
arbitrary interface location between the grid points and is not limited to middle 
point interfaces (Patankar and Spalding, 1972). An algorithm for a semi- implicit 
method (Patankar, 1980) was used as the procedure for modeling the flow field. 
The time step is determined automatically based on the size of the grid points 
and the velocity of flow fields. A detailed description of the computational 
algorithm applicable to microgravity fluid management is illustrated in our 
earlier studies (Hung et al., 1990 a,b,c). Figures 5(A) and 5(B) show the 
distribution of grid points for the dewar tank with probe for GP-B dewar 
container in the radial -axial plane and radial -circumferential plane, 
respectively, in cylindrical coordinates. 

IX. Numerical Simulations of Sloshing Dynamics 
for the Operating GP-B Spacecraft 

Computer simulations of sloshing dynamics applicable to orbital GP-B 
spacecraft have been carried out during the time period of the present research 
contract supported by the NASA grant NAG8-938. Results of the research are 
published in the following scientific journals: 

(1) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Effect of the Baffle 
on the Spacecraft Fluid Propellant Viscous Stress and Moment 
Fluctuations, Transactions of the Japan Society for Aero^ 
nautical and Space Sciences . 35, 187-207, 1993. 

(2) Hung, R. J., and Shyu, K. L. , Suction Dip, Liquid Residual and 
Slosh Wave Excitation During Liquid Draining in Microgravity, 

Advances in Space Research . 13(7) . 147-154, 1993. 

(3) Hung, R. J., and Shyu, K. L. , Liquid Settlement, Resettlement, 
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Slosh Wave Excitation and Geyser Motion During Reorientation in 
Microgravity, Advances in Space Research . 13(7) . 155-163, 1993. 

(4) Hung, R. J., Lee, C. C. , and Leslie, F. W. , Effect of the Baffle 

on the Asymmetric Gravity- Jitter Excited Slosh Waves and Spacecraft 
Moment and Angular Momentum Fluctuations, Journal of Aerospace 
Engineering . ( United Kingdom ) , 207 . 105-120, 1993. 

(5) Hung, R. J., and Pan, H. L. , Asymmetric Slosh Wave Excitation in 
Liquid-Vapor Interface Under Microgravity , Acta Mechanics Sinica . 
9(4) . 298-311, 1993. 

(6) Hung, R. J., and Lee, C. C., Effect of the Baffle on On-Orbit 
Spacecraft Fluid System Angular Momentum Fluctuations and Action 
of Viscous Stress Moment Exerted on the Container, Proceedings of 
National Science Council (A), 17, 214-235, 1993. 

(7) Hung, R. J., and Pan, H. L. , Differences in Gravity Gradient and 
Gravity Jitter-Excited Slosh Waves in Microgravity, Transactions 
of the Japan Society for Aeronautical and Space Sciences . 36, 153- 
169, 1993. 

(8) Hung, R. J., and Shyu, K. L. , Liquid Resettlement and Slosh Wave 
Excitation During Fluid Reorientation in Microgravity, Acta 
Astronautics . 27 . in press, 1994. 

(9) Hung, R. J., and Pan, H. L, , Gravity Gradient or Gravity Jitter 
Induced Viscous Stress and Moment Fluctuations in Microgravity, 
Fluid Dynamics Research . 14(1) . 29-51, 1994. 

(10) Hung, R. J., and Pan, H. L. , Cryogenic Helium System Angular 
Momentum and Moment Fluctuations Driven by Gravity Gradient 
Acceleration in Microgravity, Advances in Cryogenic Engineering . 



32, 231-240,1994. 

(11) Hung, R. J., and Lee, C. C. , Effect of the Baffle on the Cryogenic 
Helium Container Angular Momentum and Moment Fluctuations Due to 
Asymmetric Gravity Jitter Excited Slosh Waves, Advanc es in Cryogenic 
Engineering . 22, 261-270, 1994. 

(12) Hung, R. J., and Pan, H. L. , Gravity Gradient Induced Viscous Stress and 
Moment Fluctuations of Cryogenic Helium System in Microgravity, Advances in 
Cryogenic Engineer ing . 22 . 241-250, 1994. 

(13) Hung, R. J., and Shyu, K. L. , Cryogenic Liquid Hydrogen Suction Dip and 
Slosh Wave Excitation During Draining in Microgravity, Advances in Cryogenic 
Engineering . 39 . 219-229, 1994. 

(14) Hung, R. J., and Lee, C. C., Effect of the Baffle on the Cryogenic Helium 
Container Viscous Stress and Moment Fluctuations in Microgravity, Advances 
in Cryogenic Engineering . 39 . 251-260, 1994. 

(15) Hung, R. J., Lee, C. C., and Leslie, F. W. , Dynamic Characteristics of the 
Partially Filled Rotating Dewar of the Gravity Probe -B Spacecraft, Acta 
Astronautics , 32(3) . 199-209, 1994. 

(16) Hung, R. J., and Lee, C. C., Effect of Baffle on Gravity Gradient 
Acceleration Excited Slosh Waves and Associated Viscous Stress Force 
Activated Spacecraft Dynamic Fluctuations, Aeronautical Journal (Royal 
Aeronautical Society of United Kingdom), in press, 1994. 

(17) Hung, R. J., Lee, C. C., and Leslie, F. W. , Dynamics of Gravity Probe-B 
Spacecraft Due to Gravity Jitter Induced Cryogenic Helium Disturbances in 
Rotating Dewar, Advances in Space Research . 14(5) . 115-122, 1994. 

(18) Hung, R. J., Pan, H. L. , and Leslie, F. W. , Fluid System Angular Momentum 
and Moment Fluctuations Driven by Gravity Gradient or Gravity Jitter in 
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Microgravity, Journal of Flight Sciences and Space Research (Zeitschrif t 
fur Flugvissenschaf ten und Weltraumforschung ) . 18 . 195*202, 1994. 

(19) Hung, R. J., and Pan, H. L. , Mathematical Formulation of Gravity Gradient 
and Gravity Jitter Accelerations Acting on the Fluid Systems of On-Orbit 
Orbital Spacecraft, Proceedings of National Science Council (A) . 18(3) , 236- 
247, 1994. 

(20) Hung, R. J., and Pan, H. L. , Numerical Modeling of Bubble Oscillations in 
Microgravity, Scientific Computing and Automation . 10(6) . 21-25, 1994. 

(21) Hung, R. J. , and Pan, H. L. , Turn-Around Gravity Jitter Excited Slosh Waves 
and Their Effects on Angular Momentum and Viscous Stress Exerted on 
Spacecraft, Proceedings of National Science Council . (A) . 18 . 353-370, 1994. 

(22) Hung, R. J., and Shyu, K. L. , Liquid Hydrogen Suction Dip and Slosh Wave 
Excitation During Draining Under Normal and Reduced Gravity Environments, 
Transactions of the Japan Society for Aeronautical and Space Sciences . 36 . 
225-248, 1994. 

(23) Hung, R. J., Pan, H. L. , and Long, Y. T., Sloshing Dynamics Modulated 
Cryogenic Helium Fluids Driven by Gravity Gradient or Jitter Accelerations 
Associated with Slew Motion in Microgravity, Acta Mechanica Sinica . 10, in 
press, 1994. 

(24) Hung, R. J. , and Lee, C. C. , Effect of Baffle on Gravity Gradient Accelera- 
tion Excited Slosh Waves in Microgravity, Journal of Spacecraft and Rockets . 
31 . in press , 1994 . 

(25) Hung, R. J., Pan, H. L. , and Long, Y. T., Peculiar Behavior of Helium II 
Disturbances Due to Sloshing Dynamics Driven by Jitter Accelerations 
Associated With Slew Motion in Microgravity, Cryogenics . 34(8) . 641-648, 
1994. 
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(26) Hung, R. J., Long, Y. T. , and Pan, H. L. , Sloshing Dynamics Induced Angular 
Momentum Fluctuations Driven by Jitter Accelerations Associated with Slew 
Motion in Microgravity, Transactions of the Japan Society for Aeronautical 
and Space Sciences , 37. 217-233,1994. 

(27) Hung, R. J., and Pan, H. L. , Sloshing Dynamics Modulated Fluid Angular 
Momentum and Moment Fluctuations Driven by Orbital Gravity Gradient and 
Jitter Accelerations in Microgravity, Applied Scientific Research Journal 
(International Journal of the Applications of Fluid Dynamics), in press, 
1994. 

(28) Hung, R. J., and Lee, C. C., Effect of Baffles on Imbalance Bubble 
Configuration Due to Sloshing Dynamics Driven by Gravity Gradient 
Acceleration in Microgravity, Canadian Aeronautics and Space Journal . 40, 
185-202, 1994. 

(29) Hung, R. J., and Shyu, K. L. , Liquid Hydrogen Shut-Off Geyser Excitation 
Induced by Sloshing Dynamics During Draining in Microgravity. Aeronautical 
Journal (Royal Aeronautical Society of United Kingdom), in press, 1994. 

(30) Hung, R. J., and Lee, C. C., Effect of Baffles on Sloshing Modulated Fluid 
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Acceleration in Microgravity, Applied Scientific Research Journal . 
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Driven by Gravity Gradient and Jitter Accelerations, Journal of Guidance. 
Control and Dynamics . 17 . in press, 1994. 
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(33) Hung, R. J. , and Pan, H. L. , Effect of Baffle on Sloshing Modulated Torques 
Responsed to Orbital Accelerations in Microgravity, Journal of Spacecraft 
and Rockets . 31 . in press , 1994. 

(34) Hung, R. J., Tsao, Y. D. , and Spaulding, G. F. , Gravity Effect on the 
Lymphocyte Deformation Through Cell Shape Change, Proceedings of National 
Science Council . (C) . 17 . in press, 1994. 

(35) Hung, R. J., and Pan, H. L. , Characteristics of Spin-Up Cryogenic Helium 
Bubble Encirclement Around Inner Column of Dewar Container in Microgravity, 
Advances in Space Research . 14, in press, 1994. 

(36) Hung, R. J., Tsao, Y. D., and Spaulding, G. F. , Time Sequence Evolution of 
Human Cell Deformation in Micro- and Hypergravity, Advances in Space 
Research . 14, in press, 1994. 

(37) Hung, R. J., and Long, Y. T., Baffle Effect of Cryogenic Helium Sloshing 
Dynamics In Response to Lateral Impulsive Acceleration in Microgravity, 
Cryogenics . in press, 1994. 

(38) Hung, R. J. , and Pan, H. L, Rotational Speed and Wrapping of Different Size 
Cryogenic Helium Bubbles Around a Dewar Well in Microgravity, Aeronautical 
Journal (Royal Aeronautical Society of United Kingdom), in press, 1994. 

(39) Hung, R. J., and Pan, H. L. , Simulation of Bubble Spin-Up and Spin-Down in 
Microgravity, Scientific Computing and Automation . 10, in press, 1994. 

(40) Hung, R. J., and Pan, H. L. , Mathematical Modelling of Bubble Sloshing 
Dynamics for Cryogenic Liquid Helium in Orbital Spacecraft Dewar Container, 
Applied Mathematical Modelling , in press, 1994. 
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In addition to the publication in the scientific journals, we have also 

published the research results in the national and international scientific 
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(10) Hung, R. J., and Shyu, K. L. , Liquid Propellant Residual and Shut-Off in 
Spacecraft Propulsion System, Proceedings of First Pacific International 
Conference on Aerospace Science and Technology . Vol . III . 1396-1405, 1993. 
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V 



Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


GP-B spacecraft coordinate systems with azimuth angle ^ from 
spacecraft mass center to the center of the Earth. Coordinate (x" , 
y" , z") for slew motion and coordinate (x , y, z) for fluid mechanics 
computation. 

Geometry of the GP-B dewar container with the coordinate system 
perpendicular and tangential to the container wall. (A) Geometry in 
r-z plane and (B) Geometry in r-8 plane. 

Geometry of the GP-B dewar container with baffle-boards and their 
locations. (A) Geometry in r-z plane and (B) Geometry in r-0 plane. 
The GP-B module showing main elements of liquid helium dewar, probe 
and baffle-boards. 

Distribution of grid points in the (A) Radial-axial plane, and (B) 
the radial-circumferential plane of the cylindrical coordinates for 
the GP-B dewar tank. 
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